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ABSTRACT 


It  is  been  found  that  an  externally  injected  signal  can  strongly  influence  the  steady-state 
electromagnetic  mode  in  which  a  pulsed  electron  cyclotron  maser  (ECM)  operates.  A  careful 
experimental  examination  reveals  that  the  process  by  which  this  occurs  is  similar  to  the  oscillator 
phase  priming  noted  previously  in  magnetron  work.  It  is  found  here  that  the  degree  of  mode  control 
depends  on  the  injected  signal  power  level  (relative  to  the  ECM  noise  power  level),  frequency,  and  time 
of  application.  Mode  control  is  obtained  at  power  levels  nearly  four  orders  of  magnitude  below  that  of 
the  steady-state  ECM  output  jx)wer  and  over  frequency  bandwidths  several  times  that  of  the  cavity 
resonance  band.  The  optimum  input  signal  is  a  circularly  px)larized  wave,  co-rotating  with  the  electron 
cyclotron  motion.  The  experimental  results  are  compared  with  a  quasilinear  coupled-mode  theory.  The 
theory  is  used  to  provide  analytic  predictions  of  the  temporal  mode  evolution  through  third  order 
coupling.  The  comparison  with  experiment  indicates  that  higher  order  terms  are  non-negligible  in  some 
regimes  of  operation.  However,  this  model  suitably  predicts  ECM  steady  state  operation  over  much  of 
the  operating  parameter  space  and  gives  a  qualitative  understanding  of  the  role  of  the  external  signal 
and  pre-oscillation  noise  in  the  mode  selection  process.  Mode  priming  may  enable  stable  and  efficient 
single-mode  operation  of  the  ECM  and  other  overmoded  sources  of  coherent  radiation. 
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1.  INTRODUCTION 

It  has  been  previously  demonstrated  that  pure,  single-mode  operation  can  be  achieved  in  a 
overmoded  rf  oscillator  (electron  cyclotron  maser)  via  injection  of  an  external  signal.^  The  observed 
phenomenon  showed  qualitative  agreement  with  a  simple  theory  based  on  Lamb’s  quasi-linear 
treatment  of  the  multimode  problem  in  lasers.  Here  we  present  more  detailed  experimental  evidence 
of  this  effect  and  make  a  quantitative  comparison  with  the  quasi-linear  theory  appropriate  for  the 
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electron  cyclotron  maser  (ECM). 

Sources  of  coherent  radiation  typically  suffer  from  competition  between  different  spatial 
electromagnetic  modes  when  the  operating  frequency  becomes  high  enough  that  the  corresponding  free 
space  wavelength  is  small  compared  to  the  characteristic  source  dimensions.  This  competition  has  a 
variety  of  deleterious  effects  but  primary  among  them  are  losses  in  source  efficiency,  spectral 
degradation  of  the  output  radiation,  and  perhaps  undesirable  thermal  loading  or  electrical  breakdown 
in  the  source  components.  Consequently  a  variety  of  techniques  have  been  developed  to  either  prevent 
the  excitation  or  provide  suppression  of  unwanted  electromagnetic  modes.  Interferometric  methods 
involving  a  modified  cavity  structure  have  been  used  in  lasers  and  similar  slotted  or  complex  cavities 
have  been  used  to  stabilize  lower  frequency  rf  oscillators.  Another  common  technique  at  all  frequencies 
is  the  introduction  of  selective  loss  into  the  cavity  to  resistively  damp  unwanted  modes.  These 
techniques,  though  successful  in  some  limited  cases,  have  the  common  problem  of  being  rather 
inflexible,  difficult  designs  which  often  are  only  successful  over  a  narrow  region  of  operating 
parameters.  In  addition,  pure  single-mode  operation  at  the  optimum  efficiency  predicted  by  single¬ 
mode  theory  is  rarely,  if  ever,  achieved. 

On  the  other  hand,  it  has  been  known  for  some  time  that  the  steady  state  conditions  of  a 
nonlinear  oscillator  can  strongly  depend  on  the  initial  conditions  (pre-oscillation  conditions).  The  work 
on  oscillator  phase  control  in  magnetrons  by  David"*  clearly  demonstrated  this  phenomenon  in  the  field 
of  coherent  radiation  sources.  David  showed  that  by  applying  a  small  external  rf  signal  during  the 
early  stages  of  oscillation  build-up,  the  final  relative  phase  between  the  magnetron  and  external  source 
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could  be  controlled.  The  degree  to  which  control  was  maintained  depended  on  the  relative  size  of  the 
external  signal  to  the  magnetron  pre-oscillation  noise  level.  The  steady  state  magnetron  frequency  was 
unaffected  by  the  external  source,  and  hence  no  phase  locking  occurred  (evidently  the  external  power 
level  was  orders  of  magnitude  below  that  required  for  phase  locking  at  any  given  frequency  separation 
between  the  magnetron  and  the  external  source).  This  effect  was  termed  phase  priming  in  an  obvious 
analogy  to  the  priming  of  a  pump.  More  recent  experiments  have  verified  David's  results  in  the 
magnetron^'  ^  and  the  TCM.^ 

The  identifying  features  of  a  priming  phenomenon,  as  defined  by  this  early  magnetron  work, 
are  as  follows:  (1)  the  external  signal  has  an  observable  effect  on  the  steady  state  oscillation;  (2)  the 
external  signal  only  has  this  effect  if  it  is  applied  in  the  early  stages  of  oscillation  growth  (the 
oscillator  ignores  the  tiny  external  signal  during  the  steady  state,  and  hence  there  is  no  change  in  the 
operating  frequency  or  output  power);  (3)  the  effect  is  one  of  degree  (as  opposed  to  the  definite 
threshold  associated  with  the  onset  of  phase  locking);  (4)  the  degree  of  priming  depends  on  (a)  the 
power  level  of  the  external  source  relative  to  that  of  the  pre-oscillation  noise,  (b)  the  frequency 
separation  between  the  external  signal  and  oscillator  (this  dependence  is  weak  and  set  largely  by  the 
quality  factor  of  the  cavity). 

Experiments  involving  mode  selection  by  an  external  signal  were  performed  on  CO2  TEA  and 
Nd:YAG  lasers.® '  ^  The  phrase  injection  seeding  was  coined  to  describe  the  priming  of  modes.  Since,  as 
will  be  shown  here,  this  seeding  phenomenon  matches  all  four  characteristics  identified  in  the  earlier 
work,  one  might  more  aptly  use  the  term  mode  priming.  In  addition,  it  will  be  shown  that  the 
oscillator  phase  is,  in  fact,  primed  coincidentally  with  the  oscillator  mode. 

The  new  results  presented  here  are  based  on  experiments  carried  out  at  the  Naval  Research 
Laboratory  on  a  5  GHz  electron  cyclotron  maser.  The  electron  cyclotron  maser  mechanism  involves  the 
interaction  of  a  cyclotron  wave  on  an  electron  beam  with  a  co-propagating  transverse  electromagnetic 
wave.  The  negative  mass  instability  enables  efficient  emission  of  radiation.  The  experimental 
configuration  is  given  in  Fig.  1.  The  interaction  region  consists  of  a  cylindrical  waveguide  near  cutoff 
with  a  1.75  cm  radius  and  approximately  10  free  space  wavelengths  in  length.  Discontinuities 
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terminate  the  cavity  at  the  electron  gun  and  collector  ends.  The  electron  beam  parameters  are  20  kV  at 
2  A  with  a  measured  perpendicular-to-parallel  velocity  ratio,  a,  of  -  1.4  (with  an  accuracy  of  ~  20%). 
Competition  is  observed  between  axial  TEjj^  modes  at  the  fundamental  and  TEjj^  modes  at  the  second 
harmonic  of  the  cyclotron  frequency.  Here  q  is  the  axial  mode  index.  Capacitively  and  inductively 
coupled  probes  establish  the  presence  or  absence  of  modes.  Because  the  probes  are  coupled  weakly  to 
the  cavity,  internal  ohmic  losses  cause  the  quality  factors  of  the  lowest  four  axial  modes  all  to  be  ~ 
1200.  It  is  clear  that  there  will  be  significant  mode  competition  if  the  transit  time  broadening  is  larger 
than  the  inter-mode  frequency  spacing.  For  axial  modes  in  an  ECM  the  criterion  is 
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where  COj^  is  the  frequency  of  mode  n,  L  is  the  cavity  length,  and  p^  is  the  ratio  of  the  axial  electron 
velocity  to  the  speed  of  light  c.  From  this  relation  it  can  be  seen  that  since  the  modes  become  farther 
apart  in  frequency  as  the  axial  mode  index  increases,  then  the  mode  competition  decreases.  On  the 
other  hand,  the  mode  competition  increases  with  p^,  because  of  transit  time  broadening,  and  increases 
with  COj,  and  L  because  of  closer  mode  spacing.  Note  that  the  effect  of  increased  transit  time  by  increase 
of  L  (hence  decreasing  the  transit  time  broadening  and  therefore  the  degree  of  mode  competition)  is 
overcome  by  the  effect  of  L  on  the  mode  spacing.  In  this  experiment,  the  right  hand  side  of  the  above 
expression  is  ~  15,  which  is  much  larger  than  the  spacing  between  the  low  order  modes. 

The  general  philosophy  behind  these  experiments  is  as  follows.  The  electron  beam  is 
modulated  with  the  priming  signal  via  electron  cyclotron  resonance  absorption.  If  the  beam  modulation 
is  at  a  frequency  close  to  that  of  the  desired  mode,  then  that  mode  will  be  given  an  advantage  in  the 
ensuing  mode  competition  which  develops  as  the  electron  beam  excites  all  modes  within  the  cyclotron 
resonance  band.  There  is  no  a  priori  justification,  however,  for  expecting  such  a  simple  sequence  of 
events.  In  particular,  the  modes  may  not  be  strongly  coupled  or  one  mode  may  excite  rather  than 
suppress  another.  It  is  these  concerns  which  motivate  the  development  of  a  valid  theoretical  model. 
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It  is  found  experimentally  that  mode  control  can  be  obtained  in  a  variety  of  different  situations. 
When  the  unperturbed  ECM  oscillation  skips  between  two  or  more  modes  from  pulse  to  pulse,  any  of  the 
modes  can  be  selected  by  a  priming  signal  at  the  appropriate  frequency.  Also  oscillation  in  a  single 
stable  mode  can  be  switched  to  another  nearby  mode  over  a  restricted  range  of  operating  parameters. 
Finally,  a  case  is  studi  d  where  a  bistable  system  (two  simultaneously  oscillating  modes)  is  reduced  to 
a  single  mode  because  of  the  priming  signal.  Other  experiments  show  that  optimal  mode  priming  occurs 
when  the  input  signal  is  launched  as  a  circularly  polarized  wave  which  co-rotates  v.'ith  the  electron 
cyclotron  motion.  Using  the  degree  of  mode  control  as  a  measure  of  coupling  between  the  electromagnetic 
wave  and  the  electron  beam  it  is  found  that  the  ratio  of  coupling  factors  for  the  right-  and  left- 
circularly  polarized  waves  is  near  that  predicted  from  theory.^*^ 

Using  symbolic  logic  manipulation,  analytic  results  are  obtained  for  the  linear  growth  rate  and 
self-  and  cross-saturation  coefficients  for  any  circular  mode  in  the  ECM.  These  complex 

coefficients  then  are  used  in  the  rate  equations  for  the  temporal  evolution  of  amplitude  and  phase  of 
each  mode.  This  description  of  the  mode  evolution  is  displayed  in  phase  plots  (as  originally  done  by 
Lamb^  )  to  reveal  the  basins  of  attraction  of  the  various  stable  equilibria.  The  pre-oscillation  noise 
and  the  priming  signal  can  be  approximately  included  in  this  description  as  initial  conditions  from 
which  the  system  evolves.  The  case  considered  here  is  a  simple  one:  two  modes  interacting  solely 
through  their  amplitudes.  It  is  found  experimentally  that  such  a  situation  does  exist  for  particular 
values  of  the  ECM  opjerating  parameters.  This  experiment  is  compared  quantitatively  with  the  quasi- 
linear  theory. 

The  paper  is  organized  as  follows.  In  Sec.  11  we  review  the  quasi-linear  multimode  theory.  In 
this  approximation,  the  rate  equations  are  derived  and  the  coefficients  of  these  equations  are 
determined  analytically  up  to  third  order.  In  Sec.  Ill  the  experimental  evidence  is  presented  for  mode 
priming  which  includes  the  dependencies  on  power,  frequency,  polarization,  and  time  of  application  of 
the  input  signal.  Section  IV  contains  an  interpretation  of  the  experimental  results  and  includes  a 
quantitative  comparison  between  theory  and  experiment.  The  conclusions  are  given  in  Sec.  V. 
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II.  QUASI-LINEAR  THEORY 

A.  Rate  equations  for  the  electromagnetic  modes 

The  development  of  the  quasi-linear  theory  of  mode  competition  is  due  to  Lamb  and  has  been 
applied  to  the  ECM  by  Nusinovich.  Here  we  merely  state  the  assumptions  and  summarize  the 
procedure  for  including  the  external  signal  and  pre-oscillation  noise.  The  goal  is  to  arrive  at  equations 
revealing  the  dynamics  of  the  phase  and  amplitude  of  the  electric  field  of  the  various  modes.  Starting 
from  a  Slater  expansion  of  the  cavity  fields:^^ 


E  =  SanW®n  .  H=Xbn(t)hn  , 


Cn=  Vx 


lq(z)  2  ,  hn=  Vxen^ 


+  C.C.  =  5„(tl 

^  n 


+  C.C. 


a  wave  equation  for  ipfl)  can  be  obtained.  Here  and  are  the  orthonormal  eigenmodes  of  the 
cavity,  f^(z)  is  the  axial  profile  function  (the  z  axis  is  directed  along  the  cavity  axis,  in  the  direction 
of  the  electron  beam  flow),  Vy/i/  is  the  transverse  cavity  eigenfunction,  ©„  and  kp  are  the  frequency  and 
wavenumber  of  the  free-running  (single-mode)  oscillator  in  the  nth  mode  (n  is  an  abbreviation  for  the 
three  indices  m,  I  and  (f).  The  instantaneous  frequency  of  oscillation  is  ©„  +  yp.  The  rate  equations  for 
the  slowly  varying  amplitude  Ap(t)  and  phase  electric  field  of  the  nth  mode  can  be 

obtained  by  separately  equating  the  real  and  imaginary  parts  of  the  wave  equation  for  ap(t).  In  doing 
so,  it  is  assumed  that  Ap  and  ^'p  vary  little  over  a  period  of  the  field  oscillation  2n/ (Op.  Furthermore, 
the  good  cavity  assumption  1/Qo  -  1/Qe  <  <  1  is  used  where  the  Ohnaic  wall  loss  and  external  quality 
factors  are  given  by  Qo  and  Qe-  If  the  ac  part  of  the  beam  current  density  is  also  small,  then  the  rate 
equations  can  be  written 
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-^^n+- 

dt  2 


(On-^Re{L}-^^ 
Ap  CiJn 


fOr 


Re(Xn)  =  0 


^An  +  ^lm{L}  +-^^^lm{x„)  =  0 


Here  cOp  is  the  cold  cavity  resonant  frequency  of  the  nth  mode  and  L  is  defined  by 


L=I 

n' 


(1+i)kn  1 

+ 


kp'  Qq 


Q. 


nn' 


(1) 


(2) 


where  Q  and  Q  g*'  are  the  cross-coupling  quality  factors  between  modes  n  and  n'.  The  electronic 
susceptibility  Xp  is  defined  as 


Xp  —  ■  • 


:  Jn  (t) 


COpeo  An 


(3) 


where  Jp(t)  is  the  slowly  varying  amplitude  of  the  ac  electron  beam  current  density: 


Jlr.t)  =  +  c-cj  en(r) 

2  n 


(4) 


It  is  interesting  to  note  that  the  external  priming  signal  can  be  introduced  to  the  oscillation  in  two  ways. 
If  the  signal  is  coupled  to  the  electron  beam  before  it  enters  the  oscillator  (pre-modulation)  then  the 
current  density  of  Eq.  (4)  must  contain  a  component  which  oscillates  at  the  drive  frequency.  This 
method  of  injection  has  been  shown  to  be  preferential  for  ECM  phase  locking  and  phase  priming  since 
driver  isolation  can  be  maintained,  optimal  coupling  to  the  electron  beam  can  be  achieved,  and 
amplification  of  the  drive  occurs.  If  the  signal  is  coupled  directly  to  the  oscillator  cavity  then  a 
modification  must  be  made  to  the  surface  integrals  in  the  wave  equation.  For  the  purposes  of  the  current 
discussion,  the  external  signal  is  included  as  an  initial  condition  for  the  rate  equations  and  is  not 
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consistently  incorporated  into  the  equations.  This  approximation  simplifies  the  resulting  analysis  and 
is  justified  as  long  as  the  external  signal  is  so  small  that  the  oscillation  rapidly  grows  to  such  a  size 
that  it  is  unaffected  by  the  external  signal.  Since  /A^  >  >  (cOjj  is  the  priming  frequency)  then 

it  is  clear  that  many  cycles  of  the  imposed  signal  will  pass  before  the  oscillation  will  be  able  to 
outgrow  the  priming  signal.  For  this  reason  there  is  a  frequency  dependence  to  the  priming  phenomenon. 
The  theory  discussed  here  is  unable  to  predict  this  dependence. 

The  electronic  susceptibility  of  Eq.  (3)  depends  on  the  Fonder  components  of  the  current  density 
and  hence  on  the  electron  dynamics.  Thus,  to  solve  for  ,  one  must  solve  the  electron  equations  of 
motion  which,  in  turn,  depend  on  the  electromagnetic  fields.  If  the  electron  transit  time  through  the 
cavity  is  very  rapid,  then  the  field  quantities  Ap  and  'Fp  remain  approximately  constant.  In  this  way 
the  electron  dynamics  are  adiabatically  eliminated  from  the  problem  and  the  electron  current  density 
is  written  as  functions  of  Ap  and  'Fp.  The  current  density  can  be  expanded  in  powers  of  the  small 
quantity 

§  _  ^  I  ^orbitl  2  71  Tl 


where  e,  mg,  r^  ,  and  vj^^  are  the  electronic  charge,  mass,  Larmor  radius,  and  perpendicular  (to  the 

static  magnetic  field)  velocity  at  the  cavity  entrance,  respectively.  The  electric  field  magnitude  on  the 
electron  orbit  is  denoted  by  I  I  •  The  requirement  5  <  <  1  implies  that  the  electron  loses  only  a 
small  fraction  of  its  total  perpendicular  kinetic  energy  in  a  period  of  its  cyclotron  motion.  This 
requirement  is  typically  satisfied  for  a  resonant  interaction  (such  as  the  ECM)  which  requires  the 
electron  pseriod  to  be  well  defined,  thus  implying  many  electron  cycles  in  the  interaction  region. 

The  susceptibility  can  now  be  written  in  a  power  series: 

^  n  -  ^  {x  Ao  e* ‘  +  ‘I'o  -  ttno  +  ZZZ  Ao  Av  Ap  „  e^pp  + ... 

^  ^  \  O  <J  V  p 


(5) 
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where 


<D 


ovpn  ■ 


=  co„-aj..+  M  -0) 


and  the  summations  are  over  all  possible  modes.  If  the  frequency  separation  between  modes  is  much 
larger  than  the  mode  linewidths,  then  Eqs.  (1)  and  (2)  can  be  averaged  over  the  fast  time  scale  of  the 
intermode  phase  beating.  If,  in  addition,  there  are  no  degenerate  modes  and  the  modes  are  unequally 
spaced  in  frequency,  then  the  amplitude  equation  becomes  independent  of  the  mode  phases.  Expanding 
the  time-averaged  electronic  susceptibility  to  third  order  using  the  typical  notation  for  the  growth 
rate  and  saturation  coefficients: 


ttn  + 


Pn  An  +  S  Yn  v  Ay 


and  Ln  = 


(i-1) 


Qo  Qe 


ej 


An  fbn 


(6) 


where  the  replacements  Ona  ^n  /  ©nnnn  Pn  .  ^nvvn  Ynv  have  been  made.  Using  Eqs.  (6)  in 
Eqs.  (1)  there  results  a  sequence  of  rate  equations  for  the  slowly  varying  amplitudes  and  phases  of  the 
various  modes.  For  the  purpose  of  comparison  with  experiment,  it  is  converuent  to  write  these  equations 
in  terms  of  the  power  developed  in  the  cavity  rather  than  the  field  amplitude: 


^'Fn  -EWn  =  0)^  +  1  Re  [cn  -  Pn  Pn  -  XYov  ^v} 

* 

i  P„  =  Im  {a„  P„  -  Pn  P=  -  P„  SVn^,  p.) 


(7a) 


(7b) 


The  coefficients  in  Eqs.  (7)  are  related  to  those  in  Eq.  (5)  as  follows: 


an  =  ' 


(1  +  i) 


Qo  Ql 


Pn  = 


.  2  Q5 


®n  £0 


e 


nnnn  ,and 


Ynv  = 


2  Qq^ 
COy  £0 
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nwn 


(8) 
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The  Qq  in  these  formulae  originate  in  the  transformahon  from  amplitude  to  power  using  the  definition; 
Qq  =  0)  X  (stored  cavity  energy)/(power  dissipated  in  cavity  walls). 

B.  Evaluation  of  the  rate  coefficients 


The  quantities  O^nnn'  ®nvvn  found  from  the  electron  equations  of  motion 

considering  the  electromagnetic  fields  to  be  fixed.  Nusinovich  has  shown  that  the  electron  equations 
of  motion,  in  the  slow-time  scale  formulation,  can  be  written  as  an  iterative  sequence  of  equations  for 
the  perpendicular  electron  momentum  expanded  in  powers  of  6.  The  momentum  obtained  may  then  be 
used  to  find  J(r,t)  in  Eq.  (4),  which  in  turn  yields  Jn^)-  Using  Eq.  (3),  and  each  of  its  expansion 
coefficients  may  be  found.  Employing  the  approximation  that  a  cOq  L/(2  c)  >  >  1,  the  coefficients 
for  a  cylindrical  cavity  and  annular  electron  beam  are 
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Here  a  is  the  ratio  of  initial  perpendicular  to  parallel  electron  velocity,  Pj^q  is  the  raho  of  the 
perpendicular  electron  velocity  to  the  speed  of  light,  and  cOq  is  a  normalization  frequency.  The  initial 
perpendicular  momentum  is  Pj^o  electron  beam  guiding  center  radius  is  Rq.  The  eigenfunction  for 

a  circularly  polarized  TE  cavity  mode  has  been  written  ^jni  =  ^  ^  where  the 

perpendicular  wavenumber  is  the  wall  radius,  is  the  /th  zero  of  JT/ifx)]  and  the 
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constant  ~  ^  1-  Mode  ”n"  again  corresponds  to  the  cylindrical  indices  m, 

I,  and  q  which  identify  the  radial,  azimuthal,  and  axial  mode  numbers.  The  indices  of  mode  "V"  are 
M,  L,  and  Q.  The  factor  Rj,  is  defined 
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where  Iq  is  the  dc  electron  beam  current.  The  factors  labelled  with  a  "reduce"  subscript  arc  those 
calculated  analytically  using  the  symbolic  logic  manipulator  REDUCE  version  3.2.^^  The  form  of  the 
constants  has  been  given  by  Nusinovich  and  these  are  reproduced  in  the  appendix.  These  coefficients 
depend  on  the  axial  electric  field  profile  fq(z),  the  detuning  factor  Ap  =  2  (cOp  -  Cl)/ Mq),  and  the 
normalized  length  z  =  a  Pj^p  (Dp  z/(2  c).  Here  the  electron  cyclotron  frequency  is  denoted  by  Q.  These 
factors  of  Eqs.  (Al)  are  used  together  with  Eqs.  (9)  in  Eqs.  (8)  to  find  the  net  linear  growth  rate  and 
coefficients  of  self-  and  cross-  saturation  for  the  rate  equations  (7).  Note  that  the  axial  field  profile  is 
not  computed  self-consistently  in  this  approximation.  A  sinusoidal  profile  is  used  to  model  the  high 
quality  factor  modes  of  the  experiment: 


Figure  2(a)  shows  the  imaginary  part  of  Op,  which  is  proportional  to  the  linear  growth  rate, 
for  the  four  lowest  order  TE^^^^  modes.  Op  is  plotted  against  the  static  axial  magnetic  field  using 
typical  experimental  parameters  (beam  current  of  2  A,  beam  voltage  20  kV,  O.  of  1.4).  In  all 
calculations,  the  beam  perpendicular  velocity  is  assumed  constant  though  the  static  axial  magnetic 
field  is  varied.  Growth  occurs  when  the  imaginary  part  of  Op  becomes  positive.  There  are  q  regions  of 
strong  growth  for  a  TEi  ^  ^  mode.  Because  of  the  large  areas  of  overlap  (where  several  modes  grow 
linearly)  it  can  be  expect(Ki  that  there  will  be  mcxle  compjetition  in  this  device.  The  higher  order  terms 
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must  be  included  to  determine  vvhich  mode  survives  in  the  steady  state.  A  comparison  of  this  quasi- 
linear  theory  with  experiment  is  made  in  Sec-  IV. 

Figure  2(b)  shows  the  oscillator  frequency  shifts  due  to  the  real  part  of  ttp  for  the  four  lowest 
order  modes.  This  coefficient  modifies  the  frequency  of  oscillation  of  each  mode  via  Eq.  (7a).  It  can  be 
seen  that  the  frequency  separation  between  modes  is  on  the  order  of  .5  %  of  the  operating  frequency. 
Though  linear  theory  predicts  that  the  ECM  frequency  may  be  either  above  or  below  that  of  the  cold 
cavity  resonance,  in  this  experiment  it  is  found  that  the  resonant  frequency  is  always  upshifted  from 
that  of  the  cold  cavity  resonance.  It  has  been  found  that  the  higher  order  terms  do  reduce  the 
downshifting  of  the  ECM  frequency. 

C.  The  oscillator  phase  space 

Since  it  has  been  assumed  that  the  modifications  to  the  oscillator  frequency  due  to  the  electron 
beam  are  small  ('Fp  is  small),  then  the  power  rate  equations  (7b)  evolve  independently  of  those  for  the 
phases  (7a).  A  convenient  way  to  view  the  solutions  to  (7b)  is  to  eliminate  the  time  parameter  and  plot 
solutions  in  the  phase  space  formed  by  the  amplitude  (or  power)  in  each  mode.  This  was  first  done  by 
Lamb.  For  the  simple  case  of  two  modes  interacting  solely  through  their  amplitudes  he  found  three 
general  types  of  interaction.  These  he  termed  strong,  weak,  and  neutral  coupling,  respectively.  Strong 
coupling  implies  that  only  one  mode  survives  to  the  steady  state.  This  is  the  condition  most  often 
observed  in  the  ECM  with  a  closed  cavity  configuration.  Weak  or  neutral  coupling  often  involves  the 
simultaneous  oscillation  of  two  or  more  modes.  In  a  short  pulse  experiment,  it  is  not  always  clear 
whether  two  simultaneous  modes  are  weakly  coupled,  and  will  survive  indefinitely,  or  strongly 
coupled  but  take  a  long  time  to  approach  the  steady  state.  In  such  a  case  one  must  resort  to  the 
mathematical  criterion  for  strong  coupling,  which  is  ^  Pi  P2-  ^  typical  phase  diagram  for  two 

strongly  coupled  modes  is  shown  in  Fig.  3.  The  oscillator  follows  one  of  the  trajectories  shown 
depending  on  its  initial  conditions.  The  origin  is  an  unstable  equilibrium  so  that  both  modes  have 
positive  linear  growth.  Time  increases  along  each  trajectory  in  the  direction  of  the  arrows  so  that  the 
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oscillator  approaches  either  one  of  the  two  pure  mode  solutions  (P|  =  0,  P2  =  or  Pj  =  P®,^^  p2  =  0)  in 
the  steady  state.  There  is  a  separatrix  (basin  boundary)  which  separates  those  trajectories  that 
approach  mode  1  in  the  steady  state  from  those  which  approach  mode  2.  This  separatrix  plays  a 
primary  role  in  the  mode  evolution.  Though  its  equation  cannot  be  found  in  an  analytic  form,  it  does 

pass  through  two  known  points  in  the  phase  plane.  One  is  the  origin  and  the  other  is  the  intersection 
(henceforth  referred  to  as  point  1)  of  the  two  lines  Im'^aj  -  Pj  Pj  -  Pj^  }=  0,  denoted  by  the  lines  A 

and  B  in  the  figure.  As  can  be  seen  from  Eq.  (7b),  this  point  is  also  an  equilibrium  point,  but  it  can  be 
shown  to  be  unstable  for  strongly  coupled  modes. 

A  qualitative  understanding  of  the  mode  priming  phenomenon  can  be  obtained  by  using  the 
phase  plane  representation  of  the  mode  evolution.  When  the  oscillator  starts  up,  it  can  be  assumed 
that  there  is  a  broadband  noise  radiation  field  generated  predominantly  by  the  electron  cyclotron 
motion.  The  cavity  resonances  tend  to  filter  the  noise  so  that  there  remains  significant  noise  power  only 
near  the  resonance  frequency  of  each  mode.  Previous  studies  indicate  that  this  noise  is  distributed  in 

7 

amplitude  according  to  a  Rayleigh  distribution  and  is  randomly  phased. 

In  the  case  when  there  is  no  external  signal  present,  the  initial  condition  from  which  the 
oscillator  evolves  is  given  by  the  point  with  phase  plane  coordinates  (P"°'®®,  ) ,  where  is 

the  noise  power  level  in  the  ith  mode.  The  oscillation  will  proceed  to  a  pure  mode  2  or  mode  1 
oscillation  depending  on  whether  this  initial  condition  is  above  or  below  the  separatrix.  If  the  mode 
separation  is  much  less  than  the  cyclotron  emission  band  (this  is  the  usually  the  case  because  of  doppler 
broadening)  then  the  noise  power  in  two  adjacent  modes  should  be  of  similar  size.  Hence  the  initial 
condition  is  usually  near  the  45°  line  in  the  oscillator  phase  plane.  If  the  separatrix  is  also  in  this 
vicinity  then  slight  fluctuations  in  pre-osdllation  noise  may  push  the  initial  condition  back  and  forth 
across  the  separatrix  from  one  oscillator  pulse  to  the  next.  In  such  a  case  the  oscillator  would  skip  back 
and  forth  between  modes  from  pulse  to  pulse.  (Dn  the  other  hand,  if  the  separatrix  is  far  from  the  45° 
line  then  nearly  pure  mode  operation  would  result. 
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If  an  external  priming  signal  is  applied  during  the  start-up  phase  of  each  oscillator  pulse,  the 
initial  condition  from  which  the  oscillation  evolves  is  modified.  If  the  priming  signal  frequency  is 
close  to  that  of  mode  1,  for  example,  then  the  new  P|  coordinate  of  the  iiutial  condition  is  the  vector  sum 
of  the  noise  power  in  mode  1  and  the  priming  signal  power.  Since  the  P2  coordinate  of  the  initial 
condition  is  unchanged  and  the  priming  power  is  arbitrary,  the  initial  condition  can  be  varied  over  a 
wide  range  of  values.  In  particular,  the  initial  condition  can  be  pushed  across  the  separatrix,  thus 
selecting  a  new  steady  state  mode  with  only  a  small  amount  of  priming  power  (large  compared  to  the 
pre-oscillation  noise  power  level). 

The  equation  of  the  separatrix  must  be  obtained  in  order  to  predict  the  modal  evolution.  In 
practice,  since  the  noise  and  priming  power  levels  are  small  compared  to  the  steady  state  saturated 
power  level,  the  initial  conditions  of  interest  are  always  near  the  origin.  Hence  it  suffices  to  know  the 
separatrix  in  this  small  region  of  the  phase  plane.  It  can  be  shown  that  the  equation  obeys  the  power 
law  relation 

P2=C|P,)“' 

(10) 

near  the  origin.  Here  C  is  an  unknown  constant  which  depends  on  all  of  the  rate  equation  coefficients.  C 
may  be  found  by  numerically  integrating  Eq.  (7b)  backward  in  time  from  point  I  which,  as  previously 
mentioned,  must  be  on  the  separatrix.  Since  the  right-hand  side  of  Eqs.  (7b)  vanishes  at  point  I,  the 
integration  must  be  started  from  a  nearby  point  which  is  also  on  the  separatrix.  Such  a  point  can  be 
found  by  linearizing  the  rate  equations  (7b)  about  I.  Since  I  is  a  saddle  point,  the  equation  for  the 
separatrix  is  a  straight  line  in  this  vicinity.  The  starting  point  for  the  numerical  integration  is  then 
chosen  to  be  on  this  line.  A  linear  regression  of  the  numerical  data  yields  C. 

D.  Theoretical  Missed  Pulse  Fraction 
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The  degree  of  control  over  the  oscillator  modes  is  measured  by  determining  the  fraction  of 
oscillator  pulses  not  in  the  desired  mode  (missed  pulse  fraction,  F).  This  degree  of  control  can  be 
analytically  determined  if  the  mean  noise  power  and  noise  power  distribution  are  known.  Given  these, 
the  temporal  evolution  of  an  arbitrary  oscillator  pulse  is  found  by  randomly  selecting  a  power  level 
from  distributions  for  the  power  in  each  competing  mode,  then  comparing  the  resulting  initial  condition 
coordinates  to  those  of  the  separatrix  in  the  phase  plane.  After  repeating  this  process  for  a  large 
number  of  pulses,  the  fraction  F  can  be  determined. 

The  following  discussion  is  limited  to  two  competing  modes  with  a  priming  signal  applied  to 
mode  1,  the  "desired"  mode.  The  probability  density  Pp  of  the  noise  power  of  an  undriven  mode  (mode 
2)  follows  the  modified  Rayleigh  distribution: 

PF(Pd<iP2==e'N  dP2  . 

where  N  is  the  mean  noise  power  in  the  mode  (N  is  experimentally  determined).  The  driven  mode 
follows  the  distribution  Pp: 

.j  (2fPD\ 

PD|Pi)dP,  =  =  e  N  Io(-^)dP,  . 

where  D  is  the  drive  signal  power  and  Iq  is  a  modified  Bessel  function.^^  The  oscillator  will  operate  in 
the  undesired  mode  if  the  initial  condition  in  the  phase  plane  is  on  the  wrong  side  of  the  separatrix. 
For  the  typical  situation  shown  in  Fig.  3,  the  criterion  for  a  missed  pulse  is 

S(P2) 

Pl<S(P2)  or  Z=  ■  -  >  1  , 

where  Pj  =  S  (P2)  is  the  equation  of  the  separatrix.  To  find  the  probability  of  a  missed  pulse  one  needs 
the  probability  of  Z  >  1.  It  can  be  shown  that  the  probability  density  of  Z  is 
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P,(Z)dZ  =  ^ 


Pjjfx)  1 -e  N  dx  >dZ 


where  x  =  l/P^  ,  is  the  probability  density  of  x. 


{H 


P,Mdx  =  -^e-  H  mV!). 


and  S  is  the  inverse  of  S  [  i.e. ,  the  equation  of  the  separatrix  is  P2  =  S  (Pp  ].  The  missed  pulse 
fraction  is  the  probability  of  Z  >  1; 


■r 

•'z-i 


F=  Pz{Z)dZ  =  l-  Px(x)1-e  N  dx 


Near  the  origin,  the  equation  for  the  separatrix  is  given  by  Eq.  (10);  hence 


III.  EXPERIMENTAL  RESULTS 


A.  Mode  priming  the  ECM 


The  priming  experiments  involve  the  injection  of  an  external  signal  through  one  or  more  of  the 
coupling  probes  (Fig.  1)  into  the  ECM  cavity.  A  stable  sweep  oscillator  generates  the  signal  and 
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attenuators  and  phase  shifter  circuits  are  used  to  launch  waves  of  different  circular  polarization.  The 
output  from  the  ECM  generally  couples  to  all  probes  which  are  terminated  in  matched  loads.  The 
output  is  monitored  on  probes  LI  and  L2  to  detect  even  order  axial  modes  (electric  field  null  at  cavity 
midplane)  and  probes  #3  and  #4  to  detect  odd  order  modes.  The  mode  of  operation  is  determined  from 
(a)  frequency  measurements  using  a  tunable  bandpass  filter  in  the  ECM  output  line,  and  (b)  measurement 
of  the  spatial  field  profile  using  the  various  probes  along  the  cavity  axis.  More  than  40  dB  of  isolation 
is  used  to  prevent  output  power  from  the  ECM  from  reaching  the  sweep  generator.  To  further  ensure  that 
the  priming  source  is  unaffected,  a  large  drive  power  is  used  together  with  directional  couplers  and 
attenuators  to  lower  the  priming  signal  to  the  needed  size.  In  this  way  the  feedback  power  from  the 
ECM  is  small  compared  to  the  power  level  at  which  the  priming  source  operates. 

Mode  control  is  experimentally  shown  to  be  possible  under  a  variety  of  circumstances.  These 
include  (i)  switching  from  a  single  steady  mode,  (ii)  selecting  between  two  simultaneously  oscillating 
modes,  and  (iii)  selecting  between  modes  which  skip  back  and  forth  from  pulse  to  pulse.  An  example  of 
case  iii  is  given  in  Fig.  4.  The  oscilloscope  traces  show  the  signal  from  a  crystal  detector  monitoring  the 
output  power  of  the  ECM  with  beam  current  of  .77  A  and  magnetic  field  -1.87  kG.  (This  experiment  was 
carried  out  on  a  cavity  with  fewer  probes  and  a  different  electron  gun  than  previously  described.) 
When  no  priming  signal  is  applied  the  ECM  skips  between  the  TEjj j  (at  5.093  GHz)  and  TEi^3  (at  5.141 

GHz)  modes  from  pulse  to  pulse.  Throughout  a  given  pulse  there  is  no  change  in  mode.  The  exposure 
time  of  the  photograph  is  ~  2  seconds  which  corresponds  to  over  100  pulses.  Thus  both  modes  appear  in 
the  picture  at  the  same  time.  Application  of  a  small  priming  signal  at  either  the  frequency  of  the  TEj^| 
or  TEjj3  mode  allows  selection  of  the  respective  mode.  This  is  confirmed  by  frequency  measurements 
using  the  tunable  filter.  Figure  4(d)  shows  the  relative  phase  between  the  TE^^^  oscillation  and  that  of 
the  printing  signal  (priming  frequency  near  5.093  GHz).  It  can  be  seen  that  the  ECM  initially  starts  up 
in  phase  with  the  priming  signal  but  soon  wanders  off  because  of  the  slight  frequency  difference.  Thus 
the  ECM  is  phase  primed  at  the  same  time  that  it  is  mode  primed.  There  is  no  phase  locking  since  the 
priming  power  is  so  small  compared  to  the  ECM  output  power.  This  example  demonstrates  two  features 
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of  a  priming  phenomenon:  the  external  signal  has  an  effect  on  the  steady  state  oscillation  but  doesn't 
change  the  oscillation  frequency  or  p)ower  level. 

The  temporal  dependence  of  the  mode  priming  phenomenon  is  shown  in  Fig.  5.  Here  the 
external  signal  power  and  ECM  output  power  are  simultaneously  monitored.  In  the  free-running  case 
(Fig.  5(a))  the  ECM  skips  between  the  TEjj|  and  TEjj3  modes  from  pulse  to  pulse.  Both  modes  are 

apparent  in  each  picture  since  the  digitizing  oscilloscope  uses  interleaved  sampling  (samples  from 
multiple  pulses).  The  definition  of  the  trace  of  a  given  mode  is  a  qualitative  measure  of  the  fraction  of 
oscillator  pulses  in  that  mode.  Hence  in  the  free-running  case  most  of  the  ECM  pulses  are  in  the  TEj|3 

mode  because  that  trace  is  brighter.  A  50  ns  wide  priming  pulse  is  applied  to  the  ECM  about  350  ns 
before  oscillator  saturation  in  Fig.  5  (b).  The  priming  frequency  is  that  of  the  TE|j|  mode.  The  priming 
signal  slightly  increases  the  fraction  of  pulses  in  the  TE|j  j  mode  since  the  lower  trace  becomes  brighter, 
relative  to  Fig.  5(a),  and  the  upper  trace  darkens.  Moving  the  priming  pulse  about  50  ns  later  in  the 
oscillation  buildup  has  the  dramatic  effect  shown  in  Fig.  5(c).  Now  virtually  all  ECM  pulses  are  in  the 
TE^^I  mode.  Moving  the  priming  pulse  another  100  ns  closer  to  the  ECM  saturation  show  that  its  effect 

weakens  (Fig.  5(d)).  Thus  there  is  an  optimal  time  at  which  the  external  signal  should  be  applied  to 
obtain  optimal  mode  control.  If  the  priming  pulse  is  too  early  in  time,  it  will  decay  before  the  ECM 
oscillation  is  initiated.  Too  late  in  the  buildup  and  the  ECM  power  is  large  enough  to  ignore  the 
priming  signal.  This  temporal  dependence  of  mode  control  is  further  evidence  that  the  phenomenon  is 
one  of  priming.  This  feature  might  be  exploited  when  priming  a  high  power  oscillator  by  isolating  the 
priming  source  from  the  oscillator  during  the  high  power  part  of  the  oscillator  pulse. 

A  global  picture  of  mode  selection  for  the  TE||^  mode  is  shown  in  Fig.  6.  Here  the  regime  of  pure 
TEjij  mode  operation  in  both  the  free-running  and  primed  cases  are  indicated  on  a  plane  formed  by  the 
operating  parameters  beam  current  and  axial  magnetic  field(  I-B  plane).  In  the  free-running  case  a 
complicated  pattern  of  pure  modes  exists  separated  by  regions  of  pulse-to-pulse  mode  skipping. 
Application  of  the  priming  signal  allows  pure  mode  operation  if  the  priming  frequency  is  sufficiently 
close  to  one  of  the  modes.  It  can  be  seen  from  Fig.  6  that  the  regime  over  which  a  pure  TEjji  mode  can  be 

obtained  increases  by  about  40  %  when  a  priming  signal  is  applied.  The  maximum  priming  signal  power 
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used  was  less  than  20  W  (outside  the  cavity).  The  regime  of  primed  TEjjj  operation  is  found  to 
displace  free-rurming  regimes  of  stable  TEj^2'  ^114^  ^118'  ^119  operation  and  pulse-to-pulse 

mode  skipping  regimes  involving  these  and  other  modes.  Other  experiments  show  that  the  regime  of 
pure  mode  operation  of  each  of  the  higher  order  axial  modes  can  be  enhanced  by  a  priming  signal  of  the 
appropriate  frequency.  This  priming  frequency  is  very  close  to  the  frequency  the  mode  would  have,  for 
the  given  beam  current  and  magnetic  field,  if  it  could  oscillate  alone.  In  addition,  an  experiment  at  the 
second  harmonic  of  the  cyclotron  frequency,  using  TE21  q  modes  at  ~  8.5  GHz,  reveals  that  the  same 

degree  of  mode  control  can  be  obtained. 

From  these  results  it  is  clear  that  control  can  be  obtained  over  the  steady  state  mode  at  the 

boundary  between  the  excitation  of  any  two  modes  in  the  I-B  plane.  This  result  is  significant  since  it  is 

a  common  feature  of  the  ECM  that  the  regime  of  maximum  predicted  efficiency  is  near  the  boundary  of 

excitation  of  competing  modes.  It  is  to  be  emphasized  that  single-mode  operation  does  not  imply  that 

there  is  no  mode  competition.  Single-mode  operation  simply  means  that  the  mode  coupling  only  allows 

a  single  mode  in  the  steady  state.  Whether  or  not  that  mode  can  be  obtained  at  optimal  efficiency  is  an 

issue  in  which  mode  competition  plays  a  crucial  role.  Often  single-mode  theory  predicts  a  maximal 

efficiency  of  a  particular  mode  at  a  point  in  the  I-B  plane  where  a  competing  mode  is  experimentally 

found.  This  type  of  single-mode  operation  is  undesirable.  Mode  priming  may  be  a  means  of  obtaining 

the  desired  mode  in  such  a  situation.  In  addition,  it  has  been  noted  that  the  highest  efficiency  regime 

for  the  ECM  often  occurs  for  beam  currents  below  the  start  oscillation  current.  It  has  been  shown 

12 

previously  that  this  hard  excitation  regime  is  accessible  via  a  priming  signal. 

Figure  7  shows  the  minimum  drive  power  required  to  achieve  pure  mode  operation  as  a  function 
of  priming  frequency.  Here  pure  mode  operation  is  defined  to  mean  that  greater  than  99.95  %  of  all 
ECM  output  pulses  are  in  the  desired  mode.  A  magnetic  field  near  1.84  kG  and  a  22  kV,  55  A  beam  yield 
a  steady  TEj|3  mode  oscillation.  The  priming  signal  is  at  a  frequency  near  that  of  the  TEj^j  mode 

(5.074  GHz).  The  experimental  points  in  Fig.  7  generally  fall  along  an  inverted  Gaussian.  There  is  a 
slight  asymmetry  between  high  priming  frequencies  (Aco  >  10  MHz,  where  A(o  s  difference  between 
priming  and  TEj^  frequencies )  and  low  priming  frequencies  Ao)  <  10  MHz.  The  frequency  dependence  is 
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not  Lorentzian,  as  would  be  expjected  from  the  cavity  resonance  profile.  It  is  clear  that  the  change  in 
mode  control  with  drive  frequency  is  not  simply  due  to  the  finite  resonance  width  of  the  cavity.  The 
priming  power  required  for  a  given  level  of  mode  control  is  found  to  increase  exponentially  with  A(o. 
Another  unusual  feature  of  Fig.  7  is  that  the  frequency  at  which  the  least  priming  power  is  required  for 
pure  mode  operation  is  not  the  same  as  the  TEjjj  frequency.  In  this  experiment  and  in  all  other  cases 

studied  (including  experiments  at  the  second  cyclotron  harmonic),  better  mode  control  is  obtained  with 
a  priming  signal  of  frequency  slightly  above  that  of  the  free-running  oscillation.  This  same  feature  is 
present  whether  the  competing  mode  frequency  is  above  or  below  that  of  the  primed  mode.  The 
explanation  may  be  that  the  oscillation  frequency  during  the  build-up  is  not  the  same  as  that  during 
the  steady-state  (which  is  the  measured  frequency).  However,  the  theoretical  model  presented  here 
diiows  for  no  oscillator  dynamics  due  to  the  priming  signal  so  this  unusual  feature  must  be  explained 
using  a  more  detailed  description  of  the  build-up  scenario.  From  Fig.  7  it  is  seen  that,  when  priming  at 
the  optimal  frequency,  the  necessary  priming  power  (measured  outside  the  cavity)  for  pure  mode 
oscillation  is  nearly  three  orders  of  magnitude  below  the  output  oscillator  power  in  the  TEj^  mode. 

When  the  probe  coupling  is  taken  into  account,  the  priming  power  is  nearly  four  orders  of  magnitude 
below  the  oscillator  power  inside  of  the  cavity.  The  bandwidth  over  which  mode  control  can  be 
maintained  is  several  times  the  cavity  resonance  width. 

Using  a  linearly  uptapered  axial  magnetic  field  it  is  possible  to  excite  two  simultaneously 
oscillating  modes.  This  is  because  the  cyclotron  resonance  condition  varies  along  the  length  of  the 
device.  A  low  frequency  mode  oscillates  near  the  electron  gun  (upstream)  while  a  high  frequency  mode 
grows  in  the  downstream  region.  The  faster  growing  mode  is  downstream  of  the  slower  growing  one.  In 
this  way  the  downstream  mode  grows  to  a  large  size  before  encountering  any  adverse  electron  beam 
modulation  due  to  the  upstream  mode.  These  modes  are  still  strongly  coupled  but,  since  they  interact 
with  different  parts  of  the  electron  beam,  the  time  taken  to  reach  equilibrium  is  much  longer  than  the 
pulse  length.  It  is  possible  to  eliminate  the  faster  growing  mode  entirely  by  priming  the  upstream 
mode.  Priming  the  downstream  mode  has  no  effect  on  the  mode  evolution.  Reversing  the  direction  of 
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the  magnetic  field  taper  causes  the  faster  growing  mode  to  be  the  upstream  mode  and  thereby 
eliminates  the  slower  growing  mode  entirely. 

B.  Coupling  of  the  priming  signal  to  the  electron  beam 

The  degree  of  mode  control  obtained  by  priming  is  expected  to  be  related  to,  among  other  things, 
the  degree  of  coupling  of  the  input  electromagnetic  wave  to  the  electron  beam.  Thus  the  degree  of  mode 
control  can  be  used  as  an  indirect  measure  of  the  nature  of  the  electromagnetic  wave  which  participates 
in  the  electron  cyclotron  maser  mechanism.  Here  we  show  that  optimal  mode  control,  for  our  cavity, 
occurs  when  the  input  wave  is  launched  with  circular  polarization,  co-rotating  with  the  electron 
cyclotron  motion. 

The  theoretical  expression  for  the  perturbation  in  electron  motion  by  a  TE  wave  contains  terms 
proportional  to  (+  (o  -  v^  -  n  fl)  .  Here  is  the  axial  wavenumber,  v^  the  axial  electron  velocity  and 

n  the  cyclotron  harmonic  number.  The  largest  deviations  in  electron  position  are  produced  by  these 
terms  when  +  co  -  k^  v^  -  n  =  0.  Obviously  the  plus  sign  in  front  of  to  must  be  chosen  for  this 

synchronism  condition  to  be  satisfied.  This  means  that  the  electromagnetic  wave  rotates  in  the  same 
direction  as  the  beam  cyclotron  wave.  The  electrons  thus  feel  virtually  static  fields  in  the  co-rotating 
frame  of  reference  The  experimental  objective  here  is  to  show  that  regardless  of  now  the 
electromagnetic  wave  is  polarized,  upon  injection  into  the  ECM,  it  is  the  co-rotating  component  of  the 
wave  which  predominantly  couples  to  the  electron  beam.  Because  the  axial  magnetic  field  is  directed 
opposite  to  the  electron  flow,  the  best  coupling  is  expected  to  occur  via  a  left-hand  circularly  polarized 
electromagnetic  wave. 

Injecting  a  signal  (at  the  cyclotron  fundamental  frequency)  into  the  ECM  through  a  single 
capacitively  coupled  probe  excites  a  standing  TEj|  wave.  The  structure  of  this  field  in  the  cylindrical 

cavity  is  given  by 

E  =  E(j^Jl(k_Lr)  sin<|>r -h-^^^cos({)  .  (13) 
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Here  the  dependence  on  axial  coordinate  is  suppressed  for  compactness  of  presentation.  The  polar 
coordinates  are  r  and  ({)  and  the  probe  is  located  at  ({)  =  90°  (see  Fig.  8(a)).  Using  the  abbreviations 

1  clJi(kj_r) 

A  =  j^Wkxr)  a„dB  =  -^, 

1 6 

the  linearly  polarized  wave  in  (13)  can  be  written  in  terms  of  circularly  polanzed  waves: 

E  =  aJ^'I  eK“t+<!>+^/2)  ^  I  gj(cot-<t)+7C/2)  j  A  ^  ^j(cot+(l))  ^  |  gjicot-tj)) 

Here  the  right-hand  (left-hand)  component  has  a  negative  (positive)  sign  in  the  exponential.  The 
time  averaged  power  in  the  left-hand  circularly  polarized  (LHCP)  component  of  the  wave  travelling 
along  with  the  electron  beam  is  (using  the  orthogonality  of  the  right-  and  left-  hand  circularly 
polarized  waves) 

<PLH>=f^)(A2  +  B2)rd4,dr=^  ,  (15) 

s 

where  the  integration  is  over  the  cavity  cross-section.  It  can  easily  be  shown  that  the  power  into  the 
right-hand  (RHCP)  component  is  identical  to  that  of  Eq.  (15).  Thus  half  of  the  probe  input  power,  ^IN' 
goes  to  the  left-hand  wave  as  is  expected  since  the  total  wave  is  linearly  polarized. 

The  configuration  of  the  two  probe  excitation  is  shown  in  Fig.  8(b).  The  probes  are  spatially 
separated  by  90°  in  azimuthal  coordinate  and  are  excited  with  amplitudes  Ej  and  E2  and  with  a  phase 
difference  6.  The  total  drive  power  into  the  ECM  is  P|  +  P2  The  total  field  in  the  input  wave  is 

the  sum  of  Eq.  (13)  (with  E  =  e^®),  and 


E2  =  E2  (a  cos({)  r  -  B  sm<{> 


/ 


(16) 


23 


the  field  excited  by  the  second  probe.  Writing  Eq.  (16)  in  terms  of  circularly  polarized  waves,  adding  it 
to  Eq.  (14),  and  finding  the  time  averaged  power  into  the  right-  and  left-hand  waves  as  before 

\  co^  +  E2^  ±  2  E1E2  sin0]  J  (A^  -t-  B^)r  d(})  dr  ,  (17) 

s 

where  the  upper  (lower)  sign  is  taken  for  the  power  into  the  left-  (right-)  hand  circularly  polarized 
wave.  For  the  special  case  of  Ej  leading  E2  by  90°  (0  =  90°)  the  power  into  the  left  hand  wave  in 
Eq.  (17)  is 


^Plh)-2^IN  + VP1P2  (18) 

Further  specializing  to  equal  amplitudes  injected  into  the  two  probes:  (Plh)  =  I’iN '  of  i^put 

power  couples  to  the  LHCP  electromagnetic  wave.  As  long  as  the  probe  diameter  is  a  small  fraction  of 
the  cavity  circumference  a  pure  circularly  or  linearly  polarized  wave  may  be  excited.  In  the 
exjjeriment,  cold  test  results  indicate  that  the  assumption  that  a  single  probe  excites  a  pure  linearly 
polarized  wave  is  justified.  More  than  80  %  of  the  input  power  was  injected  into  the  linearly  polarized 
component. 

The  degree  of  control  of  the  modes  in  the  ECM  is  measured  by  monitoring  the  signal  from  a 
crystal  diode  morutoring  ECM  output  power.  The  ECM  is  tuned  to  a  parameter  region  where  a  single 
stable  mode  oscillates  in  the  free-running  case.  The  input  priming  signal  frequency  is  equal  to  that  of  a 
competing  mode.  As  the  drive  jX)wer  couples  to  the  electron  beam  the  ECM  begins  to  mode  skip  into  the 
driven  mode.  The  input  power  level  is  adjusted  so  that  even  when  the  drive  power  couples  optimally  to 
the  beam  the  ECM  continues  to  mode  skip  to  some  degree.  Two  methods  of  data  collection  are  used. 
First  a  simple  averaging  scheme  is  used  which  discrinrunates  the  fraction  of  missed  pulses  (those  not  in 
the  primed  mode)  on  the  basis  of  a  large  difference  in  amplitude  between  the  primed  and  free-running 


24 


modes.  The  other  method  involves  actually  acquiring  a  large  number  of  pulses  (~  1000)  on  a  digital 
oscilloscope  and  counting  the  missed  pulses  using  a  computer. 

The  averaging  method  uses  N  ECM  pulses  (N  =  256).  A  certain  fraction  of  these  pulses,  p,  will 
have  a  level  corresponding  to  the  free-running  mode.  The  remaining  fraction  of  the  pulses,  (1-  p), 
will  have  the  level  I2  of  the  driven  mode.  A  digitizing  oscilloscope  provides  the  average  A: 


A  = 


N 


pN  (l-p)N 

■  5^(1,  +  AI,|)  +  X<'2  ->■  Al2i) 
A=l  i=l 


(19) 


where  random  pulse-to-pulsc  fluctuations  in  oscillation  amplitude,  Al.,  are  included.  Equation  (19)  can 
be  rewritten 


pN  (l-p)N 

ZAI, .  Al2i 

-jqr  +  X  .  N  +  Pli  +  (l-p)l2  (20) 

i=l  i=l 

The  first  two  terms  in  Eq.  (20)  are  assumed  to  vanish.  Though  the  fluctuations  Al|  and  AI2  are  random 

Al 

from  pulse  to  pulse,  their  average  cannot  be  neglected  unless  "j”  «  1  arid  a  large  number  of  pulses  are 

included  in  the  average.  Though  the  first  condition  holds  in  the  experiment  the  second  does  not  hold  in 
the  limits  where  there  is  very  little  control  over  the  free  mode  or  the  mode  is  almost  completely 
controlled.  The  experimental  errors  are  expected  to  be  correspondingly  larger  in  these  regions.  The 
fraction  of  missed  pulses  (those  not  in  the  driven  mode)  is  therefore  p  =  (A  -  I2)  /  (Ij  ■  12)-  It  is  this 
fraction  p  which  is  the  measure  of  the  degree  of  coupling  of  the  drive  signal  to  the  electron  beam. 

Figure  9  shows  the  results  for  priming  power  injection  using  one  and  two  probes.  The  three 
experiments  shown  involve  measurement  of  missed  pulses  (using  the  averaging  method)  as  a  function  of 
priming  power.  Priming  power  is  injected  three  different  ways  and  the  results  are  shown  as  a  function  of 
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the  p>ower  into  the  LHCP  wave  in  each  case.  The  first  expjeriment  is  single  probe  drive.  The  power  into 
the  LHCP  wave  is  determined  from  Eq.  (15).  The  second  experiment  consists  of  driving  two  probes  with 
equal  amplitudes  and  90°  out  of  phase.  This  configuration  generates  a  pure  LHCP  wave.  The  last 
experiment  shown  in  the  figure  utilizes  two  probes  with  arbitrary  input  amplitudes  and  90°  out  of 
phase.  The  power  into  one  probe  is  held  constant  while  that  into  the  second  is  varied.  Equation  (18)  is 
used  in  this  case  to  determine  the  coupling  to  the  LHCP  wave. 

The  general  result  in  all  the  experiments  is  that,  as  before,  the  degree  of  mode  control  decreases 
with  decreasing  drive  power.  It  is  clear  from  Fig.  9  that  regardless  of  how  the  amplitude  of  the  drive 
signal  is  varied,  the  degree  of  mode  control  depends  predominantly  upon  the  drive  power  in  the  LHCP 
wave.  From  this,  it  can  be  inferred  that  the  electromagnetic  wave  primarily  couples  to  the  gyrating 
electron  beam,  through  electron  cyclotron  absorption,  when  the  wave  is  co-rotating  with  the  electrons. 

A  further  test  of  this  conjecture  is  provided  by  an  experiment  in  which  the  drive  power  is 
constant  but  the  phase  between  the  two  driven  probes  is  varied.  The  drive  power  input  to  each  probe  is 
approximately  equal.  The  results  of  this  investigation  are  shown  in  Fig.  10.  The  fraction  of  missed 
pulses,  p,  varies  approximately  sinusiodally  with  the  phase  angle  0.  This  is  to  be  expected  if  the 
degree  of  mode  control  is  simply  proportional  to  the  drive  power  input  in  a  LHCP  wave  since  the  power 
into  the  two  circularly  polarized  waves  in  this  experiment  is  [from  Eq.  (17)] 

<P>  =  ^  PpsI  ±  sin  0  .  (21) 


(Figure  9  shows,  however,  that  the  relationship  between  mode  control  and  LHCP  drive  power  is  not 
quite  linear.)  Of  more  importance  is  the  location  of  the  points  of  maximum  and  minimum  mode  control. 
The  first  minimum  of  Fig.  10  occurs  at  ~  130°  (probe  #3  leading  probe  #4  in  Fig.  1).  When  the  32°  line 
length  correction  is  taken  into  account  the  point  of  maximum  mode  control  is  found  to  be  ~  98°.  This  is 
quite  close  to  the  90°  interprobe  phase  shift  required  to  produce  a  pure  LHCP  wave.  The  first  maximum 
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occurs  at  a  corrected  phase  shift  of  302°,  reasonably  close  to  the  270°  required  for  excitation  of  a  pure 
RHCP  wave. 

In  conclusion,  given  a  fixed  amount  of  drive  power  the  most  efficient  means  of  coupling  to  the 
electron  beam  is  at  an  interprobe  phase  shift  such  that  a  pure  LHCP  wave  is  launched  into  the  cavity. 
Alternatively,  given  a  fixed  phase  shift  between  probes  the  amplitudes  should  be  adjusted  so  that  as 
much  drive  power  as  possible  is  in  the  form  of  a  LHCP  wave.  If  the  electron  beam  flow  is  in  the 
direction  of  the  magnetic  field  (as  is  usually  the  case)  then  optimal  coupling  should  be  obtained  with 
the  RHCP  wave. 

Theory  does  not  predict  exactly  zero  coupling  between  the  RHCP  wave  and  the  electron  beam  or 
perfect  coupling  for  the  LHCP  wave  when  the  TEj|  field  pattern  and  the  annular  electron  beam 

geometry  are  taken  into  account.  A  more  realistic  prediction  can  be  made  by  calculation  of  the  beam- 

wave  coupling  coefficient  as  defined  by  Chu.^^  The  interaction  between  the  beam  and  the 

electromagnetic  wave  is  really  only  dependent  on  the  electric  field  strength  tangential  to  the  electron 

orbit.  The  beam  power  gain  is  proportional  to  the  square  of  this  electric  field.  Expanding  the  electric 

17 

field  of  a  circular  TEjj.,j^  mode  about  an  off-axis  electron  gyro-center  yields 

E({)  ~2^ml^Tnl  Jm-l(k77i^Ro)  [  Tl)  -  12(^771/  Tl  )]  COS(CDt'(}))  , 

where  Tl  is  the  electron  Larmor  radius.  Here  only  the  first  spatial  harmonic  is  considered.  A  thin  ring 
of  electrons  with  guiding  centers  at  a  radius  Rq  is  assumed.  In  the  experiment  the  left-hand  wave  has 
eigenvalues  m=l,  n=l  and  the  coupling  coefficient  is 

^LH  =  ^l)  -  J2(^l  *■!,)]}  (22) 

The  coefficient  of  the  right-hand  wave  (m=:-l,  n=l)  is 

2 

KrH  =  {j-2(k-llRo)[l-2(k-lirL)-Io(k.ll  'l)]} 


(23) 
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Using  Rq  =  .9  cm  and  =  .23  cm  (consistent  with  the  beam  a),  the  ratio  of  the  coupling  coefficients  of 
Eqs.  (22)  and  (23)  equals  55.2.  As  expected,  the  theoretical  coupling  of  the  co-rotating  (LHCP)  wave  to 
the  beam  is  much  stronger  than  the  oppositely  rotating  wave.  A  rough  experimental  determination  of 
this  ratio  can  be  made  from  the  data  of  Figs.  9  and  10.  The  missed  pulse  fraction  can  be  eliminated  by 
transforming  p  in  Fig.  10  to  power  coupled  into  the  left-hand  wave  using  Fig.  9.  The  curve  corresponding 
to  pure  left-hand  wave  excitation  in  Fig.  9  is  used.  This  power  is  proportional  to  the  actual  power 
coupled  to  fields  on  the  electron  orbit  (Eq.  (19)  gives  the  predicted  coupling  constant).  Figure  11  shows 
the  new  plot  relating  phase  between  the  probes  to  power  actually  coupled  into  the  fields  on  the  electron 
orbit.  The  same  sinusoidal  shape  as  Fig.  10  is  apparent  here. 

The  power  on  orbit  is  produced  from  the  separate  contributions  of  the  left-  and  right-hand 
circularly  polarized  waves  (as  given  by  Eq.  (18): 

^orbit  ~  V^l  ^2  6)  +  (Pjtq  -  V^l  ^*2 

where  Kg  and  are  proportional  to  the  coupling  coefficients  of  the  right-  and  left-hand  polarized 
waves  to  the  wave  on  the  electron  orbit.  Varying  the  coefficients  Kg  and  K^  to  match  the  data  in 
Fig.  11  results  in  the  solid  line  for  a  coupling  coefficient  ratio  Kg  /  K^of  55.2.  The  curve  has  been 
adjusted  for  the  10°  systematic  error  apparent  in  Fig.  10.  There  is  some  difference  in  the  predicted 
extrema  in  the  sinusoid  but  the  overlap  with  the  experimental  points  is  fairly  good.  One  major 
uncertainty  in  this  analysis  is  that  the  averaging  method,  up»on  which  the  measurement  of  the  missed 
pulse  fraction  is  based,  suffers  its  largest  uncertainty  in  the  regions  of  very  weak  and  very  strong  mode 
control.  This  means  that  the  data  points  at  the  extrema  in  Fig.  11  are  not  very  accurate.  In  addition, 
translating  small  or  large  missed  pulse  fractions  into  equivalent  power  via  Fig.  9  is  not  precise.  For 
example,  an  uncertainty  in  the  missed  pulse  fraction  of  1  %  at  the  maximum  point  in  Fig.  10  causes 
uncertainties  in  the  coupling  ratio  of  a  factor  of  two  since  Fig.  9  gives  normalized  powers  anywhere 
between  .5  and  1.5.  It  is  found  that  variation  of  the  coupling  coefficient  ratio  by  a  factor  of  two  in 
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Eq.  (18)  still  produces  reasonable  agreement  with  the  experimental  points  in  Fig.  11.  Thus  it  can  be  said 
only  that  the  measurement  of  the  relative  coupling  strength  of  the  right-  and  left-hand  waves  is  in 
rough  agreement  with  theory. 

IV.  DISCUSSION 

A.  Mode  coupling 

The  rate  equation  coefficients  O-nre^^ce'  Pnreduce'  ^nvreduce  functions  of  three 
parameters:  cavity  length  L,  axial  wavenumber  k^,  and  detuning  between  the  field  frequency  and 

cyclotron  frequency  A.  If  the  axial  field  profile  is  assumed  to  be  that  of  the  unperturbed  cavity  mode 
(q{z),  then  k2  =  (JK/L  and  the  parameters  are  L,  q,  and  A.  Here  we  apply  these  theoretical  results  to 

the  mode  priming  experiment. 

Figure  12  shows  the  TE^n  coefficients  as  a  function  of  A  for  the  experimental  situation.  The 
subscript  on  the  coefficients  CC,  p,  and  Y indicates  the  q  index  of  the  axial  mode  (hence  TEju  indicated 
by  1,  TE^j2  f’y  etc.).  The  quality  factors  of  all  modes  are  1200,  the  electron  beam  voltage  is  19.98  kV, 
beam  current  is  .34  A,  perpendicular-to-parallel  velocity  ratio  is  1.4,  beam  guiding  center  radius  is 
.9  cm,  and  the  cavity  length  is  60  cm.  The  cold  cavity  resonant  frequencies  of  the  lowest  order  axial 
modes  are  5.033,  5.052,  and  5.081  GHz,  respectively.  The  normalization  frequency  is  5.0  GHz.  The 
detuning  parameter  Aj  in  Fig.  12  is  calculated  using  the  cold  cavity  resonant  frequency  of  the  TEm 

mode  instead  of  the  exact  field  frequency.  The  exact  frequency  is  a  function  of  the  real  parts  of  the 
coefficients  as  well  as  the  steady  state  power  levels  [Eq.  (7a)].  This  approximation  is  quite  good  and 
can  easily  be  improved  by  iteration.  Figure  12  (a)  shows  the  exact  result  for  Oj,  including  both  the 

"linear"  and  "quadratic"  bunching  terms.'’  Because  of  this,  the  growth  rate  does  not  equal  zero  when 
Aj  =  0.  For  all  coefficients,  it  can  be  seen  that  the  real  and  imaginary  parts  are  90°  "out  of  phase"; 

hence  the  real  parts  go  to  zero  near  the  points  where  the  imaginary  parts  are  at  an  extremum.  This 
implies  that,  to  this  degree  of  approximation,  the  effect  of  the  electron  beam  on  the  field  frequency  is 
minimized  when  A^  is  such  that  maximum  growth  or  saturation  occurs.  The  curves  in  Fig.  12  are 


29 


"universal"  in  the  sense  that  the  oscillator  coefficients  are  related  to  A  in  the  way  shown  even  though 
the  field  frequency  and  power  might  vary  during  an  iterative  calculation.  Figure  12  (b)  shows  pj 

calculated  using  only  the  quadratic  bunching  terms.  From  Eq.  (Al)  it  can  be  shown  that  the  imaginary 
part  of  P  must  be  an  even  function  of  A  while  the  real  part  is  odd  in  A.  In  both  Figs.  12  (b)  and  (c)  it  can 
be  seen  that  the  mode  can  be  either  suppressed  (Im(Pi)  >  0  or  Imfy|2}  >  0)  or  excited  (Im{P^}  <  0  or 
Im{Yj2)  <  0)  by  the  nonlinear  effects.  Near  resonance  (A  =  0)  the  suppression  dominates  while  for 
substantial  detuning,  A  >  .5,  there  is  growth  (Im{a|)  >  0)  and  the  nonlinear  terms  provide  further 
excitation. 

Figure  13  shows  the  coefficients  calculated  for  the  TEjj2  niode  as  a  function  of  the  detuning 
parameter  A2  (found  in  terms  of  the  cold  cavity  TE||2  resonant  frequency).  There  is  the  same  "phase" 

relationship  between  the  real  and  imaginary  parts  of  the  coefficients  as  seen  in  Fig.  12.  The  main 
difference  is  that  in  Fig.  13  there  are  two  regimes  of  TEjj2  positive  linear  growth  with  corresponding 
suppression  by  Im{P2}  and  Im{Y2i}  on  the  edges  of  the  excitation  region  closest  to  A2  =  0  and  excitation  on 
the  opposite  ends  of  the  growth  regimes.  Similar  results  are  obtained  for  higher  order  axial  modes. 

For  the  purposes  of  comparison  with  experiment,  the  imaginary  parts  of  the  coefficients  are 
shown  as  functions  of  the  magnetic  field  in  Fig.  14  for  the  lowest  three  modes.  The  region  of  positive 
linear  growth,  Im{(Xj)  >  0,  moves  to  lower  magnetic  field  for  higher  order  modes.  This  agrees 

qualitatively  with  the  experimental  observations  in  Fig.  6.  The  width  of  the  predicted  excitation 
regions  are  somewhat  different  than  that  experimentally  observed.  This  is  in  large  part  due  to  the  fact 
that  the  coefficients  are  calculated  using  the  cold  cavity  resonant  frequencies  rather  than  the  actual 
oscillation  frequencies.  All  modes  are  strongly  excited  (Im{|3)  <  0  and  Im{7}  <  0)  at  the  low  magnetic 

field  edge  of  the  first  growth  region  (1.81  - 1.83  kG).  Near  the  beam-wave  resonance  point  (-1.86  kG  for 
TEjjj,  -1.85  kG  for  TEij2'  -1-84  kG  for  TEj|3)  the  nonlinear  terms  become  positive  and  suppress  the 
oscillation.  At  higher  magnetic  fields  appear  regions  of  Im(a)  <  0  (no  linear  growth)  yet  Im{P)  and 
Im{Y)  >0  (excitation).  These  are  the  hard  excitation  regions  where  the  oscillator  is  not  linearly 

unstable  but  may  oscillate  if  an  external  "kick"  is  supplied.  Figure  14  shows  that  all  modes  have 
regions  where  hard  excitation  may  be  provided  via  TEjji  (Im{p|)  <  0  in  14  (a)  in  the  interval  1.87  - 1.89 
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kG,  Imni2i}  <  0  in  14  (b)  in  the  interval  1.855  -  1.875  kG,  Im03i)  <  0  in  14  (c)  in  the  interval  1.845  -  1.89 
kG)  or  either  of  the  other  two  modes  in  the  appropriate  regions.  This  indicates  that  mode  priming 
could  also  be  effective  in  the  hard  excitation  regions.  There  is  a  loss  of  frequency  selectivity  because  of 
the  cross-excitation  provided  by  the  Ynv  The  branches  of  TE]|^2  growth  at  high  magnetic 

field  are  generally  strongly  damped.  Figures  14  G>)  and  (c)  show  that  Im{7)  >  0  and  Im(P)  >  0  for  the 
bulk  of  those  regions  of  Im{a)  >  0  at  magnetic  fields  above  1.86  kG.  This  is  consistent  with  the 

experimental  observation  that  at  low  beam  currents  the  higher  order  modes  appear  primarily  on  the 
low  magnetic  field  side  of  the  TEj^j  excitation  region  (see  Fig.  6). 

B.  Fifth  order  quasi-linear  theory 

The  mode  priming  experiment  was  performed  at  a  magnetic  field  of  ~  1.838  kG.  For  this 
magnetic  field  many  axial  modes  have  positive  linear  growth  rates.  Table  I  shows  the  coefficients 
corresponding  to  the  experimental  parameters.  The  high  order  TE113  and  TE114  have  larger  self-  and 

cross-saturation  coefficients  than  the  low  order  modes.  It  is  found  experimentally  that  no  modes  other 
than  the  TEj^i  and  TEj-[2  can  be  detected  during  the  mode  evolution.  Injection  of  an  external  signal  at 

the  expected  frequency  of  a  weak  mode  is  found  to  generally  enhance  the  mode  amplitude.  In  this  way 
even  very  weak  modes  can  be  detected.  No  evidence  of  such  modes  was  found  in  this  experiment.  Thus  it 
can  be  inferred  that  the  high  order  modes  are  strongly  damped  and  do  not  play  a  role  in  the  resulting 
mode  evolution. 

Table  I  shows  that  all  third  order  coefficients  for  the  TEju  mode  are  negative.  This  implies 
that  all  modes  excite  the  TEm.  Since  there  is  no  saturation  of  the  TE|n  to  third  order,  it  is  clear  that 
fifth  order  terms  must  be  included.  The  terms  which  must  be  added  to  Eq.  (7  b)  are 

|PnEE8nv,PvPA  . 


-Im 
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For  the  case  of  two  competing  modes  this  correction  adds  three  new  terms  to  each  rate  equation.  The 
new  terms  in  the  TEm  power  rate  equation  are 

-Im{6iii)p?-Im(5ii2}P?P2-Im{5,22}PiP2 

Since  the  rate  equations  now  involve  second  degree  polynonnials  in  Pj  and  P2,  the  critical  points  lie 

along  second  degree  curves  (see  curves  A  and  B  in  Fig.  15)  instead  of  the  straight  lines  A  and  B  shown  in 

Fig.  3.  The  curves  still  play  the  same  role  (trajectories  have  infinite  or  zero  slope  when  they  intersect 

one  of  the  curves)  and  their  intersection  is  on  the  separatrix. 

The  stability  criteria  for  the  single  mode  oscillation  equilibria  (Pj  =  0,  P2  =  P  2  )  s^id 
(P^  =  P  ,  P2  =  0)  shown  in  Fig.  15  are 

2  2 

02  >  -Im{5222}  (Pf‘)  /  ai-Im{Yi2)Pf‘  -Im{5i22)  (pT)  <0  and 

2  2 

Ol  >  -  Im  ISiu)  (p  T  )  >  «2 '  P  f  -  Im  {8211)  (P  f  ‘  )  <  0  ^  <  25) 

respectively.  The  corresponding  criterion  for  the  instability  of  the  intersection  point  is 

Aj-j  +  A97  >  0  and  ^21  ^12  ^  ■^11  ^22  ’  (26) 

where 

=o„  -2Iml|3„)(Pr)^  -l'nll'nmlP2*-31m{6„„„)(pf  )^  -  2Im(8^)pf  P“' 

2 

-  {5nn,n,}  (p  '  n  ^  m  , 

2 

Anm  =-In'{WPf  )  -  2  Im  (5^)  P  f  P  f  ,  n^m. 


When  the  single  mode  equilibria  are  stable  and  the  intersection  equilibria  is  unstable,  then  the  modes 
are  strongly  coupled  (assuming  the  linear  growth  rate  of  both  modes  is  positive).  Inus  the  criteria  in 
Eqs.  (25)  and  (26)  are  those  for  strongly  coupled  modes  in  the  fifth  order  theoiy. 


32 


The  fifth  order  coefficients  are  not  calculated  here  but  are  estimated  from  the  experimental 
results.  The  self-saturation  is  found  from  the  experimental  single-mode  steady  state  power 

p  . 

*  n  * 


The  experimental  results  give  IrnlSm)  =  4.85  x  10^  J  and  ImfS???)  =  3.84  x  10^  J  The  two  cross¬ 
coupling  terms  6i,nm  ^nmm  Perform  simitar  roles  and,  since  neither  is  known,  arbitrarily  set 

to  zero.  The  chosen  subject  to  the  constraints  of  Eqs.  (25)  and  (26)  and  are  fit  to  the  mode 

priming  experiment  (data  of  Fig.  16).  The  values  Im{5^22)  =  1-613  x  10^  J  and  Im{82ii)  =  3-0  x  10^  ] 
are  obtained  in  this  way.  The  modal  evolution  using  the  set  of  coefficients  through  fifth  order  is  shown 
in  Fig.  15.  The  curves  A  and  B  are  ellipses  in  this  case.  The  equation  of  the  separatrix  is  found  in  the 
way  previously  described. 

C.  Comparison  of  mode  priming  theory  and  experiment 

Figure  16  shows  the  experimental  results  of  priming  a  TEjn  mode  from  a  pure  TEj  J2  naode.  The 

drive  signal  polarization  is  adjusted  for  optimal  mode  control.  As  the  drive  signal  becomes  large 
compared  to  the  mean  noise  power  level  (.3  p.W)  the  ECM  approaches  a  pure  TEi  j2  nruxle.  The  solid  line 

shows  the  results  obtained  from  the  theory  described  in  Sec.  11  D  using  the  coefficients  found  in  Sec.  FV 
B.  The  theoiy  predicts  the  missed  pulse  fraction  F  that  would  be  observed  in  a  trial  experiment 
consisting  of  an  infinite  number  of  pulses.  Since  the  actual  experimental  trial  length  is  only  256  pulses, 
a  statistical  uncertainty  enters  when  comparing  the  experiment  to  the  theory.  However  the  standard 
deviation  in  missed  pulse  fraction  of  a  large  number  of  trials  of  256  pulses  each  can  be  found  from  the 


binomial  distribution: 
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Of  = 


V 


F(1-F) 

256 


This  standard  deviation  is  indicated  in  Fig.  16  by  the  error  bars.  The  uncertainty  Op  is  small  for  low 
drive  power  because  there  is  a  large  missed  pulse  sampling  in  each  trial.  At  high  drive  power  the 
uncertainty  is  a  large  fraction  of  F,  but  is  relatively  small  (compared  to  intermediate  drive  levels)  since 
there  are  few  missed  pulses  detected  (F  is  small). 

After  taking  the  statistical  uncertainties  into  account,  the  theory  agrees  well  with  the 
experimentally  measured  missed  pulse  fraction.  This  indicates  that  the  qualitative  explanation  for 
the  mode  priming  phenomenon  set  forth  here  has  validity.  Further,  this  theory  is  adequate  for 
prediction  of  the  mode  priming  power  required  for  a  given  degree  of  oscillator  control  in  regions  of 
parameter  space  where  the  third  order  coefficients  are  positive.  A  further  improvement  can  be 
obtained  by  taking  fifth  order  coefficients  into  account.  This  will  expand  the  region  of  oscillator 
parameter  space  over  which  the  theory  has  applicability.  This  work  is  in  progress. 

V.  Conclusion 


Injection  priming  by  an  external  signal  is  shown  to  be  a  powerful  and  flexible  means  of  mode 
control  in  an  overmoded  electron  cyclotron  maser.  The  priming  technique  makes  a  much  larger  region  of 
ECM  parameter  space  accessible  to  pure  mode  operation.  The  techrdque  is  effective,  to  some  degree,  for 
all  strongly  coupled  axial  modes  which  otherwise  would  suffer  mode  skipping  or  stable  parasitic 
oscillation.  It  is  anticipated  that  this  technique  will  also  eiable  control  of  transverse  modes  as  long  as 
the  filling  factors  of  the  different  modes  are  similar. 

A  comparison  is  made  between  the  experimental  mode  priming  results  and  a  quasi-linear  theory 
of  mode  competition.  It  is  found  that  the  theory  gives  a  good  qualitative  picture  of  the  mechanism  by 
which  the  injected  signal  selects  modes  in  the  strongly  coupled  case.  It  is  found  that  for  a  broadly 
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applicable  quantitative  prediction  of  mode  control,  terms  higher  than  third  order  must  be  included  in 
the  theory. 

An  indirect  measurement  is  made  of  the  relative  coupling  of  the  right-  and  left-hand  circularly 
polarized  waves  to  the  electron  beam.  It  is  found,  as  predicted  from  theory,  that  the  wave  with  co¬ 
rotating  polarization  is  the  dominant  electromagnetic  wave  involved  in  the  electron  cyclotron  resonance 
mechanism.  The  experimental  relative  coupling  strength  agrees  with  the  theoretical  value  to  within 
experimental  uncertainties. 
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Appendix:  Rale  Equation  Coefficients 


The  coefficients  of  linear  growth  and  self-  and  cross-saturation  are 


«n™duc.=  I  W:(Z)|U,<Z)-  iVn{^}dZ  , 


Pn  -  I  w^(z)  I  2  >^*^2  )jf  dz"  6z  +  U*iz )  v^z’)  -  u^z')  viz)  (Sz 


dzjdz 


Ynv  =  -  2i 

reduGB 


1  {  1  ^  ^1  fo  ^  ^  ^ 


f 


Unlz”)  uJX)  V*v{z')  dz 


dz /dz 


where  Wn(z)  =  f<2{z)  e' ^nZ  ,  Un(z)=  (  Wn{z')dz  ,andVn{^=  (  Un(z’)  dz  . 
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TABLE  I 

Calculated  rate  equation  coefficients  corresponding  to  the  experiment:  beam  voltage  19.98  kV, 
beam  current  .34  A,  a  =  1.4,  magnetic  field  of  1.838  kC. 


TEiii  TEi12  TEi13  TEi14 


«!  = 

:  .321 

02  =  -476 

as  = .349 

04  =  .31 

Pl  = 

-.0048 

P2  =  .0092 

P3  =  -017 

P4  =  .018 

Yi2  = 

-.0088 

Y21  = .0079 

Y31  =  .0195 

Y41  =  .027 

Yi3  = 

-.0075 

Y23  = .0213 

Y32  =  .0324 

Y42  =  .0388 

7i4  = 

-.0008 

Y24  =  .0208 

Y34  =  .037 

Y43  =  -042 

HGURE  CAFIIONS 


1:  Configuration  of  the  ECM  with  inset  showing  the  midplane  probe  geometry. 

2:  (a)  Magnetic  field  dependence  of  the  ImfCX^)  for  various  modes.  Positive  ImlOj^)  means  linear 
growth,  (b)  Linear  predictions  of  the  oscillation  frequency  [from  Re  { CCj^ }  ]  as  a  function  of 
magnetic  field. 

3:  Mode  evolution  for  two  strongly  coupled  modes.  Arrows  indicate  the  direction  of  increasing  time. 

4:  Mode  priming  of  the  ECM.  (a)  Free  oscillation  skipping  between  the  TEj^  (upper  trace)  and  the 
TE||3  modes,  (b)  Pure  TEm  mode  primed  operation  (injected  signal  frequency  equals  TEm 
frequency  of  5.093  GH2).  (c)  Pure  TEj^3  mode  primed  operation  (injected  signal  frequency  is  5.141 
GHz),  (d)  Relative  phase  between  the  TEjj|  oscillation  and  priming  signal  corresponding  to  the 
experiment  shown  in  (b). 

5:  Temporal  dependence  of  mode  priming.  The  top  trace  in  each  shows  detector  measurement  of  rf 

drive  p>ower  while  the  bottom  measures  ECM  output  power,  (a)  No  drive  signal,  mode  skipping 
between  TE^^^i  (lower)  and  TEj|3  (upper),  (b)  Drive  signal  applied  too  early  in  buildup,  (c) 

Optimal  application  time  of  drive  pulse,  (d)  Drive  pulse  applied  too  late. 

6:  Enlargement  of  the  region  of  pure  TEj  mode  oscillation  due  to  a  20-W  drive  signal  at  the 
oscillation  frequency  (shaded).  Regions  of  displaced  pure  mode  operation  are  indicated. 

7:  Minimum  drive  power  required  to  prime  a  pure  TEm  mode  from  pure  TEn2  as  a  function  of 
frequency  separation  between  the  drive  signal  and  TEju  oscillation. 

8:  Geometry  for  (a)  one-probe  and  (b)  two-probe  excitation. 
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9:  Mode  control  as  a  function  of  normalized  power  into  the  left-hand  circularly  polarized 

electromagnetic  wave.  The  three  experiments  involve  driving  one  probe,  two  probes  with  equal 
power  and  90°  phase  lead  on  probe  #1,  and  two  probes  at  90°  with  unequal  power. 

10:  Degree  of  mode  control  as  a  function  of  phase  shift  between  probes.  Phase  shifts  corresponding  to 
pure  right-hand  (dot)  and  left-hand  (dash)  waves  are  indicated. 

11:  Effective  power  on  the  electron  orbit  as  a  function  of  relative  phase  between  probes.  The  line 
shows  theory  for  a  left-  to  right-hand  coupling  ratio  of  55.2. 

12:  Universal  curves  of  TEm  coefficients  for  different  detuning  parameters  ;  (a)  CX|  ,  (b)  (5| ,  (c)  7^2 
with  real  part  (dash)  and  imaginary  part  (solid). 

13:  Universal  curves  of  TEjj2  coefficients  for  different  detuning  parameters:  (a)  (X2  ,  (b)  P2'  ^21 

with  real  part  (dash)  and  imaginary  part  (solid). 

14:  Imaginary  part  of  the  coefficients  for  the  three  lowest  order  axial  modes  vs.  magnetic  field;  (a) 
TEm  ,  (b)  TE^|2  '  ^^113-  respective  scale  and  dimensions  are  a ,  (ns'^  );  P  x  2  x  10'®  ,  (]  ‘^  )  ; 

7x2x10-®  ,(J‘^). 

15:  Mode  evolution  for  the  experimental  case  including  fifth  order  coupling  terms. 

16:  Degree  of  mode  control  achieved  with  a  given  drive  power  and  .3  jiW  noise  power.  Error  bars  show 
the  standard  deviation  in  predicted  control. 
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